Experimental observation of speckle instability in Kerr random media 
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In a disordered medium with Kerr-type nonlinearity, the transmitted speckle pattern was predicted to become 
unstable, as a result of the positive feedback between the intensity fluctuations and the nonlinear dependence of 
the local refractive index. We present the first experimental evidence of speckle instability of light transversally 
scattered in a liquid crystal cell. A two-dimensional controlled disorder is imprinted in the medium through a 
suitable illumination of a photoconductive wall of the cell, whereas the nonlinear response is obtained through 
optical reorientation of the liquid crystal molecules. Speckle patterns with oscillating intensity are observed 
above a critical threshold, which is found to depend on the scattering mean free path, as predicted by theory. 

PACS numbers: 42.25.Dd, 42.55.Zz 



In a disordered medium, the transmitted intensity measured 
at a point results from the interference of light incoming from 
all possible scattering trajectories inside the medium. These 
constructive and destructive interferences give rise to inten- 
sity fluctuations [1]. The resulting speckle pattern observed 
in transmission is a fingerprint of the medium traversed by 
the wave and is highly sensitive for instance to any displace- 
ment of scatterers within the medium. This property has been 
widely used in applications such as speckle imaging [ 2] or dif- 
fuse wave spectroscopy Short and long range speckle 
correlations, which can also results from frequency or polar- 
ization shift of the incident beam, have raised considerable 
interest in the last two decades |5[|. But changes in speckle 
pattern can also be achieved by intensity variation of the in- 
put beam when considering a nonlinear scattering medium 
with a Kerr-type nonlinearity. Local intensity fluctuations 
induce local refractive-index changes, which in turn modify 
the optical trajectories, resulting in a redistribution of the in- 
tensity fluctuations within the sample and in a new transmit- 
ted speckle pattern. Under certain conditions, this feedback 
mechanism was predicted to become unstable leading to the 
spontaneous oscillations of the speckle pattern It was 
pointed out in |7|] that speckle instabilities differ from com- 
mon nonlinear instable regimes observed in nonlinear optics, 
since here the feedback is provided by the scattering. Accord- 
ing to Skipetrov and Maynard 10, 0], the instability threshold 
p = (ri2l} 2 (L / £) 3 > 1 is reached for any value of the nonlin- 
ear correction to the refractive index r^i, as long as the mean 
free path t is small enough. To date, experimental observation 
of speckle instability has never been reported. 

In this letter, we report the first experimental evidence of 
speckle instability. A liquid crystal (LC) cell serves simulta- 
neously as the scattering and Kerr nonlinear medium. A two 
dimensional (2D) computer-generated image of random scat- 
terers is projected on a photoconductive wall of the cell, gen- 
erating a controlled 2D random distribution of the refractive 
index within the LC layer. A probe laser beam is transversally 
scattered while inducing a nonlinear (NL) Kerr-reorientation 
of the LC molecules. Speckle pattern spontaneous oscillations 
are observed. Instability threshold is identified and found to 



be dependent on scattering mean free path, as predicted by 
theory. A frequency cascade is observed at higher probe inten- 
sity. Just above threshold, the oscillation frequency is found 
to be of the same order of the inverse of the NL relaxation 
time, which is interpreted in terms of NL energy transfer and 
natural gain selection. Light scattering in random and nonlin- 
ear media has attracted considerable attention recently with, 
for instance, the exploration of the impact of a nonlinearity 
on Anderson Localization Jsl [To|] or the role of disorder on 
solitons 111 in . Our experimental study is among the very first 
Jl2ll to explore the dynamic interplay between nonlinearity 
and multiple scattering. 




FIG. 1. (color online). Experimental configuration. The disorder 
configuration is imprinted within the LC (bottom-left inset) by illu- 
minating (Iw) the photoconductive wall (PH) of the LC light- valve 
with a random distribution of scatterers (top-right inset), whereas the 
nonlinearity is provided by the molecular reorientation of the liquid 
crystal under the action of the probe probe (ip). (G=glass cover). 

The experimental setup is schematically represented in 
Fig[TJ The 20mm x 30mm liquid crystal light- valve (LCLV) is 
composed of a L = 55/im-thick LC film in the nematic phase 
and a 1 mm-thick Bii2SiO20 (BSO) photoconductive sub- 
strate, sandwiched between two transparent indium-tin-oxide 
electrodes A random pattern (here a random collection 
of dark disks in a white background, as shown in top-right 
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inset of Fig. Q]) is computer-generated and transmitted to a 
36. 9mmx 27.6mm spatial-light-modulator (SLM), spatial res- 
olution 1024x768 pixels, which serves as a 2D random mask. 
The local voltage applied across the LC film decreases where 
light impinges on the BSO crystal, inducing local LC reori- 
entation. A given light intensity distribution is therefore con- 
verted into a corresponding distribution of molecular orien- 
tation. As a result of LC birefringence, a 2D refractive index 
distribution invariant along the longitudinal direction, n(x,y), 
is formed within the LC film, which reproduces the initial ran- 
dom pattern. Because the BSO crystal is preferentially sen- 
sitive to the blue-green region of the spectrum, a CW diode- 
pumped solid state laser at 532 nm is used as the writing beam 
with an intensity of a few mW/cm 2 (2mW/cm 2 ). The refrac- 
tive index distribution is monitored by a 8 -bit CCD-camera 
and is shown in bottom-left inset of FigQ] The geometri- 
cal characteristics of the LC scattering medium are fully con- 
trolled by adjusting the filling fraction, (/>, the diameter D and 
positions of the cylindrical scatterers. By shrinking the image 
on the computer screen, the actual scatterer diameter can be 
varied from D = 100/im to 1 /im, leaving (j> constant. The re- 
fractive index of the scatterers, n s , varies between n Q = 1.52 
and n e = 1.75, according to the 8-bit coded gray level of the 
SLM mask, where n Q and n e are, respectively, the ordinary 
and the extraordinary index of the LC. In the following, we 
choose (j) =40%, D = 50 or 20 /im and n s = n e = 1.75 while 
for the LC matrix n = n . 

A linearly polarized HeNe laser (632. 8nm wavelength, 
5mW power, beam- waist i^o=380 /im, intensity 1.1 W/cm 2 ) 
serves as the probe beam. Its polarization is parallel to the 
LC nematic director. The probe beam is transversally scat- 
tered as it propagates along the LC layer and is back-reflected 
at the LC/BSO interface with a total reflectivity of 22%. The 
back-reflected near-field speckle pattern is imaged onto a sec- 
ond CCD camera and is shown in Fig. [2 We checked that 
polarization is conserved. This confirms that propagation re- 
mains in the paraxial approximation limit and allows to esti- 
mate the mean free path in the transverse direction. Indeed in 
this limit, the wave equation reduces to a Schrodinger equa- 
tion where time is replaced by longitudinal coordinate, z, to 
describe a particle moving in a 2D potential —V(x,y) = 
— k 2 [n 2 (x, y) — ng]: 
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where A(x, y, z) exp(— ikn^z) is the wave field and no 
is the effective refractive index defined by no = 
S~ x j s n 2 (x, y)dxdy iTull . When considering transverse scat- 
tering, the relevant wavewavector is no longer k = uj/c but 
rather the transverse wave- vector k± = uo/v = 2/wo, with 
v/c = k/k± ~ 2 in our case. Because this ratio can be larger 
than 1, the product k±£ <C k£, making transverse Anderson 
localization of light possible in finite-sized random systems 
even for small index contrast OH Q • In our system, a 
crude estimate of the mean free path from Mie theory (infi- 
nite cylinders with index n e in a matrix n ) and single scat- 



tering approximation yields £ w 100/im (resp. £ « 40/im) 
at 632. 8nm for D = 50/im (resp. D = 20/im). The lo- 
calization length in 2D is given by £ w £ex.-p(7rk±£/2), 
which yields £ = 228,6/mi (resp. £ = 55.7/im). As it 
propagates freely along the longitudinal direction, the beam 
will extend transversally, experiencing successively a ballistic 
regime (t < To = £/v), a subdiffusive multiple- scattering 
regime (to < t < Ti oc ) Il5ll and a crossover to localiza- 
tion (t > Tioc). To identify the regime reached within the 
round-trip ballistic time % = 2L/c spent in the LC cell, 
we estimate h/ro = (d/£)(k/k±)) w 2.2 (resp. 5.5), and 
t h / Tioc = 2(t b /r )exp(-7rk ± £) w 0.8% (resp. 5.7). Al- 
though based on a rough estimate of the mean frfee path, these 
calculations demonstrate that the multiple scattering regime is 
well established, even though light may not be yet fully trans- 
versely localized. 
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FIG. 2. (color online). Temporal modulation of speckle intensity and 
corresponding power spectrum within an area encompassing a single 
speckle spot (delineated by the white square shown in the speckle 
pattern CCD image). Time-averaged intensity in this area is (a) 1=40 
bits and (b) 1=140 bits. 
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FIG. 3. (color online). Amplitude of oscillations vs average intensity 
at the different frequencies seen in Fig [5Jd, obtained by averaging 
spectra measured at positions of the speckle pattern with same time- 
averaged intensity. 

During its round-trip propagation in the cell, the 
extraordinary-polarized probe beam experiences nonlinear 
scattering as its electric field changes molecule orientation. 
Surprisingly, the speckle pattern is not static but we observe 
with a naked-eye a temporal modulation of each speckle spot 
(as shown in the movie fujll ). These oscillations are not ob- 
served for ordinary-polarized light and therefore cannot be 
attributed to thermal fluctuations. They are found to aver- 
age out when spatial integration is performed, which shows 
that no phase relation exists from spot to spot. It is impor- 
tant to point out that our observations are fundamentally dif- 
ferent from standard self-phase modulation (SPM) described 
in the literature for non-disordered media [17]. The uncorre- 
cted oscillations for different speckle spots confirm the non- 
deterministic nature of the speckle instability, in contrast for 
instance with the phase-coherent rings observed in SPM. 

We analyze the oscillating speckle behavior by monitoring 
over time the modulated intensity, spatially averaged over a 
small area including a single speckle spot. The temporal oscil- 
lations and the corresponding amplitude spectrum are shown 
in Fig. [21 for two speckle spots with different time averaged 
intensity. Low intensity speckle spot shows a single peak at 
j o=2. 39 Hz (Fig. [2^). This has to be compared to the free 
relaxation time of the reorientational nonlinearity of the LC, 
tnl = lL 2 /7r 2 K, where 7 is the LC rotational viscosity 
and K the elastic constant. For typical values of 7 = 0.02 
Pa.sandK = 15 pN El El], 

we find with good agreement 
fo w 1/tnl- The same time constant, approximately 400 ms, 
has been independently found by measuring directly the re- 
sponse time of the optical reorientational effect. The speckle 
instability is therefore solely driven by the slow response time 
of the reorientational Kerr nonlinearity. This is another reason 
why it cannot be related to thermal instabilities since thermal 
fluctuations are expected to be faster than molecular orienta- 



tion of the LC (of the order of 100 /is EH])- Besides, given 
that the intensity used is relatively small and given the small 
absorption coefficient of pure LC , we can safely neglect 
thermal heating due to laser illumination. 

Higher intensity speckle spots show several spectral peaks 
at fo/2, /o/4 and also at 2/o (Fig. [J]}). Increasing the probe 
beam intensity results in the same spectral cascade. Remark- 
ably, the intensity distribution of the speckle pattern offers a 
direct exploration of the intensity dependence of the insta- 
bility at fixed input beam intensity. This is further analyzed 
by averaging spectra measured at positions with same time- 
averaged intensity. Four peaks with different threshold are 
identified and their amplitudes relative to the spectral back- 
ground are plotted versus speckle spot intensity in Fig. [51 
Different slopes associated with different threshold are ob- 
served for each oscillation frequency. Next we investigate 
the instability threshold dependence with disorder. According 
to the theoretical predictions ELS], the instability threshold 
p = (ri2l) 2 (L/£) 3 > 1 is expected to decrease with decreas- 
ing mean free path. We compare the intensity dependence of 
the peak amplitude at the main oscillation frequency fo for 
scatterers with diameters D = 50 /im and 20 /im, correspond- 
ing respectively to £ w 100/im and £ w 40/im. Changing 
the scatterer dimensions does not affect the oscillation fre- 
quency fo. However the instability threshold is clearly seen 
in Fig. [H to decrease with increasing scattering (decreasing 
£). This clearly demonstrates the essential role of disorder in 
speckle instability and the threshold dependence on the mean 
free path, as predicted by theory. A more systematic study 
of the threshold dependence with disorder and a comparison 
with theoretical predictions are under investigation. 
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FIG. 4. Amplitude of oscillations at fo=239 Hz vs. average inten- 
sity for two different scatterer diameters D = 50 /im and 20 /im, 
corresponding respectively to scattering media with mean free path 
£ « 100/im and £ ~ 40/im. The threshold, which is materialized by 
the crossing of the lines, decreases with decreasing £. 

We interpret these experimental observations in terms of 
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parametric amplification. In a parametric process, a pump 
beam at f p and a probe beam at f m interfere to create a grating 
within the nonlinear material. The pump beam is itself scat- 
tered by this grating and nonlinearity-driven energy transfer is 
achieved from the pump to the probe when Af = f p — f m ^ 
0. This transfer of energy via the Kerr nonlinearity has been 
interpreted as an amplification process, where the gain reads 
as G(Af) = 2Af r NL [l + (Af tnl) 2 }' 1 [20]. The cor- 
responding gain curve is peaked at Af = 1/tnl- Starting 
from noise, energy from the pump at f p will be preferentially 
transferred into photons at f m = f p — 1/tnl resulting in 
the beating of the two beams. In the same way, the spon- 
taneous speckle oscillations observed in our experiment re- 
sult from the beating between the incident HeNe beam and 
the nonlinearly-generated beam at the maximum of the gain 
curve, giving rise to a temporal modulation at fo « 1/tnl- 

In conclusion, we have observed speckle instabilities in 
transverse scattering of light in random LC films with reori- 
entational Kerr effect. The instability is solely attributed to 
the combination of scattering and Kerr nonlinearity, with un- 
corrected oscillations for different speckle spots. It depends 
on polarization, excluding any thermal effect. The inten- 
sity distribution of the speckle allows a direct probing of the 
intensity-dependence of the instability spectrum. A threshold 
is found which depends on the degree of disorder (£). This 
confirms earlier theoretical predictions Q]. The oscillation 
frequency just above threshold is equal to the inverse of the 
NL relaxation time and is interpreted in terms of the beating 
between the incident beam and a generated light emitted at a 
frequency equal to the maximum of the gain curve associated 
with the Kerr effect. Two regimes of speckle instability have 
been predicted theoretically 12111 : a slow NL regime where 
tnl > td driven by the relaxation time of the non linear- 
ity, as in the experiment presented here, and a fast NL regime 
when tnl < td- This second regime was out of reach in the 
present experiment since the reorientational Kerr NL is ex- 
tremely slow. However, we conjecture that a frequency of os- 
cillation different from 1/tnl should be observed, which will 
strongly depend this time on the characteristics of the scatter- 
ing medium. As in a random laser 11221 12311 . the gain curve will 
select an eigenmode of the passive system, which in turn will 
beat with the pump. Different solutions can be suggested to 
reduce the relaxation time of LC NL reorientation, using for 
instance ferroelectric liquid crystal (r < 1ms). This would 
allow the exploration of the fast NL regime (tnl < Td), 
where the instability oscillation should depend on the scat- 



tering characteristics of the material. 
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